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ABSTRACT 

We exploit the vastly increased sensitivity of the Expanded Very Large Array 
(EVLA) to study the radio continuum and polarization properties of the post- 
starburst, dwarf irregular galaxy IC 10 at 6 cm, at a linear resolution of ~ 50 pc. 
We find close agreement between radio continuum and Ha emission, from the 
brightest H II regions to the weaker emission in the disk. A quantitative analysis 
shows a strictly linear correlation, where the thermal component contributes 50% 
to the total radio emission, the remainder being due to a non-thermal compo- 
nent with a surprisingly steep radio spectral index of between —0.7 and —1.0 
suggesting substantial radiation losses of the cosmic-ray electrons. We confirm 
and clearly resolve polarized emission at the 10 — 20% level associated with a 
non-thermal superbubble, where the ordered magnetic field is possibly enhanced 
due to the compression of the expanding bubble. A fraction of the cosmic-ray 
electrons has likely escaped because the measured radio emission is a factor of 3 
lower than what is suggested by the Ha inferred SFR. 

Subject headings: galaxies: dwarf — galaxies: individual(IC 10) — galaxies: 
magnetic fields — radio continuum: galaxies 



Introduction 



Dwarf galaxies are the most abundant type of galaxy in the local Universe ( IMarzke fc da Costa 



1997 ) and are the closest analogues to the building blocks fo r galaxy formation at high red- 
shift in ACDM m odels (see, for example. lLagos et al. 2009 ). Furthermore, in the "down- 
sizing" paradigm ( )Hu et al.lll995l : iHomeier et al.ll2006l ) dwarf galaxies increasingly dominate 
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the star^forming universe over time ( Sawicki fc Thompson 20061 ) . with s tar formation (SF ) 
perhaps being particularly enhanced within dwarfs over the past 2 Gyr ( jOrban et al.l 120081 ) . 
Stellar winds and supernova explosions of the most massive stars can easily lead to outflows 
due to their weak gravitational potential. This implies that not only does the metallicity of 
the interstellar medium increase with cosmic time, but additional ly the intergalactic medium 
becomes enriched at an early stage ( jOppenheimer fc Pavel 120061 ). Moreover, dwarf galaxies 
are supposed to have played a fundamental role in amplifying magneti c fields in the earl y 
Universe through outflowing plasma escaping from star-forming regions ( IChyzy et al.ll201ll ). 



Cosmic rays play an important role in galactic outflows, as they consist of relativistic 
particles that do not quickly lose their energy, like the hot gas. They spiral along magnetic 
field lines emitting highly linearly polarized, non-thermal (synchrotron) emission providing 
a useful tool to study magnetic field structure. This magnetic field influences how cosmic 
rays propagate. For example, nearby edge-on galaxies with radio haloes show a character- 
istic "pinched at the waist" structure in polarised emission (and so metimes even in total 
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2009; 


Soida et al. 


2011) 



Resolved radio continuum st udies of dwarf galaxies have been lacking thus far, because 
the objects are intrinsically faint. IChyzy et al.l ( 120031 ) took advantage of the excellent sensi- 
tivity of the 100-m Effelsberg telescope to study IC 10 in the radio continuum, both in total 
intensity and polarimetry. They found extended, diffuse emission coincidin g with Ha emis- 
sion, a nd polarized emission close to the non-thermal bubble described by lYang fc Skillman 
( 1 19931 ). However, they were limited by their coarse resolution (If 1) that prevents any de- 
tailed comparison with ancillar y data. A higher reso lution preliminary VLA map at 6 cm 
was published bv Ictivzvl fcoosh . Ictivzv et all fcoooh . using the VLA, found in NGC4449 
despite its weak rotation a surprisingly regular magnetic field showi ng hints of spiral struc- 
ture and "fans" possibly connected to outflows. iKepley et al.l ( 120101 ). again using the VLA, 
studied the starburst dwarf NGC 1569 which has a strong radial field, likely shaped by the 
outflowing gas. 

With the new capabilities of the expanded very large array (EVLA) and much improved 
sensitivity we can study the weak extended emission in dwarf galaxies. In this letter we 
present first results from observations of IC 10 with the EVLA at C-band, using the full 
2 GHz of bandwidth that was available at the time of observation. Additional observations 
at L-, X-, and Ka-band, along with a full interpretation of the data will be published in a 
forthcoming paper. 



IC 10 is a perfect candidate for an in depth study of a dwarf galaxy at radio continuum 
wavelengths because it is relatively bright and nearby, being a member of the Local Group 
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( jSanna et al.ll2008l ). It has a high star formation rate and a high density of Wolf-Rayet 



stars and is regarded to be in a starburst phase (IMassey fc Holmesl 12002 : ICrowther et al. 



20031 ) . The aftermath of massive star formation leaves its imprint on the ISM in the form 
of expanding bubbles and shells that can be traced in H I and Ha. IWilcots fc Miller! ( 119981 ) 
found holes in H I and Ha emission wit h large filaments in Ha extending up to a few 
hundred parsec away from the main body faunter et allbnod l. suggesting outflows. The H I 
emission is asymmetric with respect to the major axis. A feature visible in velocity maps 
suggests ongoing accretion of neutral hydrogen in the outer parts of the galaxy. There is 
a detection of diffuse hot X-r ay emitting gas in th i s gala xy that is aligned with the main 
stellar bar dWang et all l2005h . lYang fc Skillmanl (Il993h found a region with high radio 
continuum brightness and id entified it as a non-thermal bubble prior to the break out. 
Lozinskaya fc Moiseevl (120071 ) speculated that it could be generated by a hypernova rather 
than by a collection of supernovae. 

As IC 10 is close to the Galactic midplane (b = — 3?3), estimates of its distance are 
uncertain due to the unknown level of absorption. Here, we use a conservative distance of 
1 Mpc in lin e with the most recent results that i ndicate a likely distance at best slightly less 
than IMpc dKim et al.lhood : ISanna et al.ll2008h . 



2. Observations 

We observed IC 10 with the NRAO0 EVLA under the resident shared risk program 
(RSRO; project ID: AH1006). Our observations were taken in D-configuration on August 
16, 2010. We used 2 IFs each consisting of 8 sub-bands with a bandwidth of 128 MHz each. 
The sub-bands contain 64 channels at 2 MHz resolution resulting in a total bandwidth of 
2 GHz. We placed our IFs between 4.5—5.4 GHz and 6.9—7.8 GHz, respectively. We observed 
for 4h to ensure good (u, ^)-coverage and sensitivity to extended emission. Every 15min 
we observed a phase cal ibrator (J01 2+58 24); we used 3C48 as our flux density (primary) 
calibrator, adopting the iBaars et al.l ( 119771 ) flux scale. In order to measure the polarization 
leakage (D-terms) we observed the unpolarized source 3C84. The polarization angle was 
calibrated using 3C48. 

We calibrated the data with the Common Astronomy Software Applicationspackage 
(CASA) as outlined in the continuum tutorial of 3C391 on the CASA homepage!. In all 
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subbands, leakages of 5% — 15% before calibration were reduced to 0.2% (variation across 
channels) after calibration. R-L polarization angles varied by 5° — 10° with a constant slope 
across channels before calibration, and went down to 1° — 2° variations across channels after 
calibra tion. We constructed Stokes L Q ) and U images using the Multi-Scale Clean algo- 
rithm (ICornwelll 12008k iRich et al.ll2008l ) in CAS A that models extended emission accurately, 
but assumes a flat sky-spectrum. In Stokes I the integrated flux density was 131 mJy, with a 
peak brightness of 11.3 mJy. The integrated polarized intensity is 3.0 mJy with peak bright- 
ness of —58 /xJy and 47 /xJy in Stokes Q and [/, respectively. The target rms noise calculated 
for this observation with ^3 hours on-source was 4 /iJy. The rms noise in the Stokes Q and 
U images was 5 /xJy, but the Stokes I image was affected by deconvolution errors at the 
40 /iJy level. These errors were at the expected level for continuum imaging when ignoring 
spectral structure (average dynamic range ^300 for a spectral index of —0.7 across a 50% 
bandwidth, 4.5 - 7.8 GHz). 

Stokes I imagi ng was repeated with the Multi -Scale Multi-Frequency Synthesis (MS- 
MFS) algorithm ( jRaull2010l : iRau fc Cornwellll201ll ) which simultaneously solves for spatial 
and spectral structure during wide-band image reconstruction, and reached an rms noise 
of 15 /xJy. This wide-band model was then used to perform one round of amplitude and 
phase self-calibration, and the resulting image (shown in greyscale in Fig. [I]) achieved the 
rms noise of 5 /xJy only 20% higher than the theoretically expected level. We note that the 
rms is identical for Stokes /, Q ) and U ) so that we are not dominated by the signal-to-noise 
ratio in Stokes /; the self-calibration is sufficiently accurate. 

A radio spectral index image was produced by MS-MFS as well, and a post-deconvolution 
wide-band primary beam correction was applied to remove the effect of the frequency- 
dependent primary beam. One of the goals of this RSRO project was to test the accuracy 
of this approach. Therefore, this spectral index map was compared with a low-resolution 
spectral index map constructed within AIPS from primary-beam corrected Stokes / im- 
ages at each sub-band and smoothed to the resolution of the lowest sub-band. The results 
agreed, confirming that for this source, wide-band imaging via MS-MFS followed by a post- 
deconvolution primary beam correction was able to recover the same spectral information as 
traditional methods, but at a much improved angular resolution. 

The Stokes Q and U images were taken into AIPsE, and maps of the polarized intensity 
were constructed using comb with option pole that applies a correction to the noise bias. We 
also created the maps of the polarization angle in comb using the option pola. We rotate 
the polarization (E-)vectors by 90° to show the orientation of the magnetic field, where we 
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neglect any contribution of Faraday rotation. We tested whether the structure of Stokes Q 
and U changed with frequency by making images for the 16 individual sub-bands and found 
the change to be small, indicating negligible Faraday rotation within our frequency range. 

We calculated the degree of polarization and found values between 5 — 20%. These values 
are at least a factor of 5 larger than w hat we would expec t from instrumental polarization 
which is around 1% for the VLA beam ( ICondon et al.lll998l ). This was further checked by us 
as we carried out separate test observations with 3C84 at the half-power point in C-band 
and found values there of at most 3%. However, because of beam smearing and the fact that 
our source lies well within the half-power point, we expect instrumental polarization to be 
considerably smaller. We therefore conclude that our measured polarization is dominated 
by the intrinsic polarization of our source. 



Results 



Figured] shows our total power continuum radio map at 6 cm with a resolution of 9" 4 x 
7! 3. Contours start at 3 x the rms noise level and increase by a factor of about 2 with 
every further contour. We choose a logarithmic grey scale to show the noise floor that is 
evenly distributed over the full extent of the map. To our knowledge this is the deepest 
radio continuum map of a dwarf galaxy yet m ade, going deeper than the preliminary VLA 
map, also at 6 cm, published by IChyzyl ( 120051 ) which r eaches a noise of 9 u Jy (Chyzy, priv. 
comm.). The deep VLA map of NGC 1569 at 6cm by lKepley et al.l (120101 ) has an rms noise 
of 8.5 /xJy beam -1 but at a smaller beam size of ~ 4" and is much less sensitive to extended 
emission. Although there are no artefacts visible in the figure that could be ascribed to 
missing short spacings, it can not be ruled out that we miss low-level, extended emission. We 
find a total flux of 131 mJy whereas we expect, based on interpolating single dish observations 



at 10.45 GHz and 2.64 GHz (IChvzv et al 
to 30% in total power. 



2003. 



201 ll ) 190 mJy. We could thus be missing up 



We compare in Fig.[2]the radio continuum emission with a deep Ha map of lHunter fc Elmegreen 



(120041 ) . The Ha is shown as a logarithmic grey-scale map to enhance weak emission features 
away from the main stellar body. The radio and Ha emission follow each other closely. The 
radio continuum peaks at the location of the Ha emission maxima that trace H II complexes. 
Moreover, thanks to the exquisite dynamic range, we can compare the radio emission with 
extended Ha emission, including various shells and filaments. The large filament in the West 
with a length of 600 pc is the most prominent feature but we are sensitive enough to also 
detect the much weaker filaments in the Eastern half. The most prominent hole in the radio 
continuum is found at RA = 00 h 20 m 23 s , DEC = 59°18f4 and is fully contained, both in Ha 
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and radio. It was identified by IWilcots fc Millerl ( 119981 ) to be the most prominent hole in 



H I emission with hints of expansion. 

The close correlation between the radio continuum and Ha emission suggests a large 
fraction of thermal radio emission from free-free emission. This is borne out by the radio 
spectral index map presented in Fig. [31 The thermal emission has a radio spectral index of 
a = —0.1 (S oc z/*), whereas the non-thermal synchrotron emission has a spectral index of 
—0.7 (provided the cosmic ray electrons are young, otherwise the spectrum is yet steeper). 
We find that almost all compact resolved sources within the galaxy have a flat spectral 
index between —0.1 and —0.2 indicative of a dominant thermal component. These are the 
compact H II regions with intense star formation visible in Ha emission. We note that our 
radio spectral index map obtained via MS-MFS has a better resolution than that obtained 
with the conventional method, where the resolution is limited by the low frequency map. 
This allows us to reliably identify the compact sources in the radio spectral index map that 
otherwise would be confused in the lower resolution version. 

Away from the H II regions the radio spectral index steepens to values between —0.7 
and —1. This can only be explained with a non-thermal synchrotron component. Also, 
this requires synchrotron and inverse Compton losses to be important for the cosmic-ray 
electrons as the radio spectral index is significantly steeper than —0.7, the value expected 
for young cosmic-ray electrons that have been accelerated in supernova shock waves. This 
result is surprising as cosmic rays in starburst dwarf galaxies are supposed to be advected 
in the outflow and their residence time would be short compared to their cooling time. 
The radio spectral index steepens even more dramaticall y at the south-eastern e dge (RA 



= 00 h 20 m 30 s , DEC = 59°16f5) where the superbubble of lYang fc Skillmanl (119931 ) resides. 
There the values are between —1 and —1.2 indicating dominating synchrotron and inverse 
Compton losses. 

The map of the structure of the regular magnetic field is presented in Fig.[H We show the 
orientation of the large-scale magnetic field as measured from the linear polarization, where 
the length of the vectors is proportional to the polarized intensity. The degree of polarization 
in this region is 10 — 20% indicating a significant ordered magnetic field component possibly 
enhanced by compression by the expanding superbubble. The magnetic field orientation 
is approximately NE-SW, where the polarized emission extends to the southern tip of the 
galaxy. 

We also detect a magnetic field at the location of the two brightest H II regions in 
the southeastern part. At this position, away from the superbubble, the magnetic field 
orientation is more aligned with the stellar body. We can not decide, based on the current 
data, if the magnetic field traverses the H II regions, because the field could also lie in front 
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of or behind them, in a thick gas layer as is expected for dwarf galaxies. At this resolution 
we do not detect a magnetic field that extends the full extent of the galaxy. We note that 
the location of our polarized intensity and the orientation of the magnetic fi eld agrees well 
wi th that meas ured with the 100-m Effelsberg telescope (IChyzy et al.l 120031 ) and the map 
of IChyzyl ( 120051 ). lending confidence to the performance of the EVLA and the CAS A data 
reduction package. 



4. Discussion 



We compared the flux densities of the radio emission and Ha emission averaged in 
50 pc (10") boxes. The radio and Ha maps were convolved with a Gaussian to a HPBW 
of 50 pc prior to averaging. The result is present ed in Fig. [51 The Ha map was corrected 
for foreground absorption using E(B — V) = 0.75 (IBurstein fc Heilesl 11984 ). We clipped the 
radio data below 5 x the rms noise level, to ensure we only use reliable data points and 
to moderate the influence of missing flux density on our interferometric observations. A 
power law fit in the double logarithmic plot results in a power law index of 0.97 ±0.03 and is 
therefore consistent with a strictly linear relation between the two variables. The dashed line 
shows the relation if the radio flux density were p urely thermal free-free emission calculated 
from the Ha flux density (e.g. iDeeg et al.lll997l . Eq. 3), assuming an electron temperature 
of 10 4 K. The thermal flux density alone is a factor 2 lower than the measured radio flux 
density. This is consistent with a thermal fraction of 50% in this galaxy as indicated by the 
flat radio spectral index. 

The linear relation between the radio and the Ha brightness means that we can use 
both equally well as s tar form a tion t racers. The star formation rate (SFR) from the radio 
continuum is given by ICondonl (119921 ) and can be converted to a star formation r ate density. 
If we convert the Ha emission into a star formation rate by the relation given in iKennicutt 
(119981 ). we can compare the two measurements. The radio star formation rate is a factor of 3 
lower than that implied by Ha. This means that IC 10 lies beneath the radio-SFR correlation 
and is radio dim. Further studies are required to investigate how this corresponds to the 
thermal and non-thermal fraction of the radio emission. In theory we expect that the high 
frequency observations are a better tracer for the thermal component, especially since they 
require no correction for internal or foreground extinction. One goal of this project is to 
establish how well we can use radio emission as star formation tracer even in such extreme 
cases as starburst dwarf galaxies. 



Leroy et al.l (120051 ) found that the molecular depletion time is short in IC 10 and that 



the inferred SFRs are too high for the molecular gas content. This may be explained by a 
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starburst althou gh it is not clear how much gas is hidden a nd not traced by the CO sub-mm 
rotational lines ( IMadden et al.l Il997l : iBolatto et al.l 120001 ) . As in many dwarf galaxies the 
knowledge of the molecular gas content is still sparse, the development of an alternative 
SFR tracer becomes even more important than for normal spiral galaxies. 

The magnetic field structure in this galaxy lacks a strong large-scale component that 
would be apparent in radio continuum polarization. This fits the expectation that the 
standard a^-dynamo ne eds shear to o perate as near-solid body rotation will provide only 
weak amplification (e.g. IWidrowl 120021 ) . Such conditions favour studies of local magnetic 
field amplification mechanisms that are normally difficult to observe in spiral galaxies, as in 
these systems shear-enhanced disk magnetic fields dominate. The non-thermal superbubble 
in IC 10 offers the possibility to study what role cosmic rays and magnetic fields play in 
the disk-halo interface in porous interste l lar m edium conditions. State of the art MHD 
simulations by Ide Avillez fc Breitschwerdtl ( 120051 ) suggest that the magnetic field gets swept 
up in the walls of the shells but can not prevent a superbubble break-out. Eventually the 
full set of observations will allow us to get a handle on the energetics of the superbubble and 
put these models to the test. 
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Fig. 1. — Total power radio continuum at 6 cm at 9'.'4 x 7" 3 resolution as a quasi-logarithmic 
grey-scale image. Contours are at 3, 6, 10, 20, 40, 80, 150, 300, and 600 x 5/xJybeam -1 . 
The rms noise level is at 5 /xJy beam -1 . The total power map is corrected for primary beam 
attenuation. 



- 12 - 




00 20 45 30 15 00 

RIGHT ASCENSION (J2000) 



Fig. 2. — Same as Fig. [T] but as an overlay on Ha shown as a quasi-logarithmic grey-scale 
image. 
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Fig. 3. — Radio spectral index map shown as a color-coded image at 9"4 x 7'!3 resolution. 
The contours show the radio continuum emission and are similar to those in Fig.[TJ The 
frequency dependence of the primary beam was removed in a post-imaging step. 
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Fig. 4. — Polarized intensity and magnetic field orientation overlaid on the fractional polar- 
ization (grey scale) at 15" resolution. Polarized intensity contours are at 3, 6, 10, and 20 
x 7/uJy beam -1 . The vector length is proportional to the polarized intensity, where 10" is 
equivalent to 22/iJybeam -1 (corrected for primary beam attenuation). 
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Fig. 5. — Comparison of the flux densities of the radio continuum emission at 6 cm with that 
of Ra at a scale of 50 pc. The solid line shows a fit to the data as described in the text. 
The dashed line shows the expected contribution based on Ha to the radio continuum via 
thermal free- free emission (see text for details). The horizontal line indicates 5 x the rms 
noise level in the radio map, where we clipped the radio data. 



